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Abstract. We present a preliminary phylogeny for 19 species of African Agama lizards based on a maximum likelihood 
phylogenetic analysis of 1,181 bp of mitochondrial DNA sequence data. Monophyletic radiations of species in East, South 
and West Africa are supported, as well as a clade containing two species (4. doriae and A. sankaranica) distributed across 
the Sahel region. West African populations of A. agama are paraphyletic with respect to A. finchi from westernmost Kenya, 
providing further evidence for a biogeographic corridor between West and East Africa. Populations of A. agama form 
four phylogeographic groups, which suggests that 4. agama may be composed of multiple independent evolutionary lin- 


eages. 
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INTRODUCTION 


African lizards in the genus Agama are among the most 
diverse and widespread terrestrial squamates in Africa, 
making them an ideal group for testing biogeographic hy- 
potheses and conducting comparative ecological and evo- 
lutionary studies. However, the robust phylogenetic 
framework for Agama that is needed to meet these goals 
is currently lacking. Constructing a phylogeny for Aga- 
ma requires detailed investigations of phylogenetic rela- 
tionships across multiple spatial and temporal scales. In 
general, the higher-level relationships within the Agamı- 
dae are poorly understood (MACEY et al. 2006), and dense 
sampling within Agama and other African and Asian taxa 
is necessary for testing the monophyly of and identifying 
the sister taxon to Agama. At the population level, sever- 
al widespread and polytypic species, such as Agama aga- 
ma in West Africa and A. lionotus in East Africa, have un- 
clear species limits (BOHME et al. 2005) and may be com- 
posed of multiple independent evolutionary lineages. 


These broadly distributed taxa pose some difficult species 
delimitation problems (WAGNER 2007; WAGNER et al. 
2008a, WAGNER et al. 2008b) and require detailed phylo- 
geographic investigations. 


The goals of our collaborative research on Agama lizards 
are to combine our sampling cfforts from separate biodi- 
versity fieldwork programs in different countries across 
Africa to 1) infer a comprehensive molecular phylogeny 
for the genus Agama, 2) resolve species limits in wide- 
spread, non-exclusive taxa, 3) conduct phylogeographic 
investigations of the A. agama and A. lionotus complex- 
es, and 4) describe new specics, where warranted. In this 
paper, we present a preliminary phylogeny for the genus 
Agama based on an analysis of mitochondrial DNA se- 
quence data (mtDNA). In addition, we investigate the phy- 
logeography of A. agama across West Africa. 
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MATERIAL & METHODS 


A total of 19 species of Agama were included in the phy- 
logenetic analyses (Table 1). We includcd 15 specimens 
of A. agama from 10 countries across West Africa to in- 
vestigate the phylogeographic relationships among pop- 
ulations. We selected four species from the Agaminae to 
represent outgroup taxa, including Acanthocercus atricol- 
lis, Landakia stellio, Trapelus mutabilis, and Xenagama 
taylori. We rooted our phylogeny with Landakia stellio, 
although we note that the higher-level relationships with- 
in the Agaminae are poorly understood (MACEY et al. 
2000; 2006). 


DNA was extracted from tissues following the PureGene 
Animal Tissue DNA Isolation Protocol (Gentra Systems, 
Inc). Two portions of the mitochondrial genome were PCR 
amplified and sequenced. including a small fragment of 
the 16S rRNA gene (16S) and a portion of the ND4 pro- 
tein coding gene (ND4) and the adjacent histidine, serine. 
and leucine tRNA genes. Oligonucleotide primers used for 
PCR and sequencing are provided in Table 2. The 165 
rRNA gene was PCR amplified for 30 cycles (95°C 30s, 
58°C 30s, 72°C 50s) and the ND4 was amplified for 35 
cycles (95°C 30s, 55°C 40s, 72°C T min). PCR products 
were purified using ExoSAP-IT (USB Corp.). Cycle se- 
quencing products were ethanol precipitated, and then se- 
quenccd using an ABI 3730 automated sequencer. 
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Fig. 1. Phylogenetic relationships within the genus Agama based on a maximum likelihood analysis of mtDNA sequence data 
using the GTR+T model of nucleotide substitution. Maximum likelihood bootstrap values 250% are shown. 
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Table 1. Voucher numbers and loeality data for specimens used in the study. All sequenees are deposited in GenBank (Accessi- 
on Nos. GU128430 — GU128502). Voucher abbreviations are as follows: CAS, California Academy of Scicnces; MVZ, Muscum 
of Vertebrate Zoology; ZFMK, Zoologisehes Forsehungsinstitut und Museum; MCZ, Museum of Comparative Zoology; LSUMZ, 
Louisiana State University Museum of Natural Scienee; ULM, University of Louisiana at Monroe Muscum of Natural History; 
AMNH, American Museum of Natural History; AMB, Aaron M. Bauer personal eolleetion; specimen numbers beginning or en- 


ding with “X”, “T°, or “TR”, personal collector numbers for authors from the Muséum national d’Histoire naturclle. 


Species 

Agama aculeata 
Agama agama 
Agama agama 
Agania agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama agama 
Agama anchietae 
Agama armata 
Agama atra 
Agama boulengeri 
Agama caudospinosa 
Agama doriae 
Agama finchi 
Agama hispida 
Agama impalearis 
Agama kaimosae 
Agama lionotus 
Agama mwanzae 
Agama planiceps 
Agama rueppelli 
Agama sankaranica 
Agama spinosa 
Agama weidholzi 


Acanthocercus atricollis 


Laudakia stellio 
Trapelns mutabilis 
Xenagania taylori 


Locality 


Botswana, Kgalagadi District 
Cameroon 

Cameroon, Mende, monts Takamanda 
Cameroon, Yaounde 

Cameroon, Yaounde 


Equatorial Guinea, Bioko Island, Malabo 


Gabon, Barrage de Tehimbcle 
Ghana, Greater Acera Region, Tesano 
Ghana, Northern Region, Buipe 
Ghana, Volta Region 

Guinea, Diaragbela, Niger River 
Mahi 

Mauritania, Sehbabı 

Niger, Niamey 

Nigeria, Cross River National Park 
Tchad, Bol 

Namibia, Khorixas District 
Tanzania 

South Afriea, Northern Cape Provinee 
Mauritania, Nouakehott Dist. 
Kenya, Naru Moru 

Camcroon, Darc Ville 

Kenya, Malaba 

South Africa, Northern Cape Province 
Mauritania, Nouakehott Distriet 
Kenya, Nandi 

Kenya, Nakuru 

Kenya, Mara 

Namibia, Outjo Dist., Kamanjab 
Somaha, Awdal Region 

Ghana, Volta Region 

Djibouti, Dikhil 

Senegal 

Uganda, Rukungin District 
Turkey, Antalya Province 

Egypt 

Somalia, Galbed Region 


Voucher 


MVZ 198076 
ZFMK 83761 
3853 

MCZ 184561 
MCZ 184562 
CAS 207953 
ZFMK 73185 


LSUMZ H20336 
LSUMZ H20085 


MVZ 249617 
ULM 200 


AMNH 109799 


ZFMK 76838 
MVZ 23889] 
MVZ 253099 
2906 I 

AMB 7582 
ZFMK 84990 
CAS 193435 
MVZ 235764 
ZFMK 83662 
4218 X 
ZFMK 83653 
AMB 4800 
MVZ 235766 
ZFMK 83660 
ZFMK 83646 
ZFMK 82076 
AMB 7638 
MVZ 241340 
MVZ 249656 
MVZ 236458 
TR 481 

CAS 201726 
MVZ 230213 
ZFMK 64395 
MVZ 241356 


Contiguous DNA sequences were aligned and edited us- 
ing Sequencher v4.8, and multiple sequence alignments 
were generated using Muscle v3.6 (EDGAR 2004). The 
open reading frame for the ND4 gene was identified us- 
ing MacClade v4.08 (MADDISON & MADDISON 2005). The 
16S rRNA alignment included tndel-rich loop regions that 
could not be aligned unambiguously and were excluded 
from the phylogenetic analysis. Phylogenetic relationships 
were inferred using maximum likelihood using RAxML 
v7.0.4 (STAMATAKIS 2008), and implemented the GTR+T 
model of nucleotide substitution. Support values for ın- 
ferred relationships were estimated from 100 non-paramet- 
ric bootstrap replicates. 


RESULTS 


The maximum likelihood analysis of the combined mtD- 
NA sequence data (1,181 characters) supports the mono- 
phyly of the genus Agama (bootstrap = 100%; Fig. 1). The 
outgroup taxa Xenagama tavlori and Acanthocercus atri- 
collis form a strongly supported clade (bootstrap = 100%) 
that is sister to Trapelus mutabilis. Within Agama, a clade 
containing A. impalearis and A. spinosa is sister to all re- 
maining Agama, although this relationship is not accom- 
panied by strong support (bootstrap <50%; Fig. 1). Some 
regional species assemblages are monophyletic, including 
a clade of Southern African species (A. acnleata, A. an- 
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Fig. 2. Phylogeographic structure of Agama agama across sub-Saharan Africa. Samples of Agama agama are labeled by coun- 
try in the unrooted phylogeny. The maximum likelihood phylogenetic analysis places Agama planiceps sister to Agama agama. 


chietae, A. atra, A. armata and A. hispida) and a clade of 
East African species (A. kaimosae, A. mwanzae, A. liono- 
tus, A. candospinosa and A. rueppelli), each of which are 
strongly supported (bootstrap = 99% and 100%, respec- 
tively; Fig. 1). Strong support is also provided for a clade 
of species distributed across the Sahel region of sub-Sa- 
haran Africa, including A. doriae and A. sankaranica 
(bootstrap = 100%). Agama planiceps is sister to a clade 
containing all 15 specimens of 4. agama and A. finchi, al- 
though this relationship is weak (bootstrap < 50%; Fig. 
1). The interrelationships among these major clades are 
not supported by bootstrap values > 50%, except for a 
strongly supported clade (bootstrap = 100%) containing 
the Sahel species (A. doriae and A. sankaranica) and A. 
planiceps, A. agama, and A. finchi (Fig. 1). 


Populations of Agama agama in West Africa are para- 
phyletic with respect to A. finchi from westernmost Kenya 
(Fig. 1). Populations of A. agama form four phylogeo- 


graphic groups (Fig. 2). The geographic proximity be- 
tween samples is a poor predictor of phylogenetic rela- 
tionships, and three of the four clades overlap spatially 
(Fig. 2). One of the most well sampled phylogeographic 
groups within A. agama extends across West Africa from 
Mauritania to Gabon (Fig. 2). Two other clades with broad 
distributions include a group extending from Guinea to 
Cameroon and another extending from Tchad to Kenya (A. 


finch?) (Fig. 2). Finally, a clade with a more restricted dis- 


tribution includes populations from southeastern Nigeria, 
western Camcroon, and Bioko Island (Equatorial Guinea) 
(Fig. 2). 


DISCUSSION 
This study represents the first detailed molecular phylo- 


genetic investigation of the genus Agama. Monophyletic 
radiations of Agama lizards occur in East Africa, South 
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Table 2. Primers used to amplıfy and sequence mitochondrial DNA. 


Gene Primers 7 5 | Reference 
16S 16Sf: CGCCTGTTTAACAAAAACAT This study 
16sR: CCGGTCTGAACTCAGATCACGT 
NDA ND4:CACCTATGACTACCAAAAGCTCATGTAGAAGC AREVALO et al. (1994) 


LEU: ACCACGTTTAGGTTCATTTTCATTAC 


Africa, the Sahel region, and even West Africa (Fig. 1), 
but reconstructing the biogeographic history of the group 
is impeded by several factors. First, the inter-relationships 
among the major lineages of Agama inferred by the mtD- 
NA data are not accompanied by strong support, which 
would make any biogeographic scenarios based on the cur- 
rent phylogeny speculative. Second, our taxonomic cov- 
erage 1s not comprehensive and lacks approximately 10 
species that are distributed throughout Africa. Finally, the 
sister taxon of Agama remains unclear. Despite these chal- 
lenges, there are strong biogeographic signals in the cur- 
rent phylogeny (Fig. 1). For instance, the phylogenetic 
placement of A. finchi from western Kenya within the 
West African A. agama complex provides evidence for a 
biogeographic corridor between West and East Africa, a 
result that provides further evidence for a close biogeo- 
graphic affinity between these regions (WAGNER et al. 
2008c). In addition, the close relationship between A. plan- 
iceps from Namibia and the West African A. agama com- 
plex provides support for a biogeographic connection be- 
tween West and South Africa (Figs 1 and 2). This relation- 
ship is not surprising given the close similarity in mor- 
phology and colour pattern shared between these species 
(BOHME et al. 2005). 


The intraspecific phylogeny for Agama agama highlights 
the taxonomic problems that this widespread taxon intro- 
duces to Agama lizard systematics. This species 1s para- 
phyletic with respect to A. finchi (Figs 1, 2) and appears 
to be composed of multiple distinct clades. The argument 
could be made that A. finch? should be synonymized with 
A. agama to retain a monophyletic A. agama, but A. finchi 
is clearly distinct based on morphology and coloration. 
The major phylogeographic groups found within 4. aga- 
ma based on mtDNA may represent distinct species, and 
represent a good starting point for testing species bound- 
aries with multiple nuclear markers. Determining if gene 
flow is absent among these mtDNA groups is an impor- 
tant next step in resolving the systematics of the A. aga- 
ma complex. The paraphyly of A. agama also underscores 
the need for additional detailed comparative morpholog- 
ical investigations that do not rely solely on adult male 
coloration characteristics (see GRANDISON 1968). 


The phylogeographic relationships within the Agama aga- 
ma complex exhibit a surprising amount of spatial over- 
lap between clades. Phylogeographic studies of other Aga- 
ma lizards, including A. impalearis (BROWN et al. 2002) 
and A. atra (MATTHEE & FLEMMING 2002) recovered more 
typical phylogeographic patterns whereby populations 
formed geographically exclusive clades. The spatial 
overlap of mtDNA clades ın the A. agama complex could 
be evidence of the presence of multiple distinct species, 
which co-occur throughout West Africa and are as of yet 
morphologically cryptic. Conversely, this spatial pattern 
could reflect the recent expansion of distinct lineages that 
were formerly restricted to exclusive geographic areas. 
Distinguishing among these different hypotheses of pop- 
ulation history await the addition of more specimens from 
throughout West A frica and the collection of nuclear DNA 
data. This work is currently underway. 


The collaborative nature of this research project is pro- 
viding an opportunity to conduct detailed investigations 
of the phylogeny and phylogeography of Agama lizards 
than would have otherwise been impossible. Our current 
work benefits from having expanded taxonomic sampling 
and multiple individuals of each specics, more detailed 
outgroup sampling, dense population sampling in A. aga- 
ma and A. lionotus species, and the addition of nuclear 
DNA sequence data. 
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